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ABSTRACT 

We investigate structural properties of dark matter halos of disk galaxies in hierarchical ACDM cosmology, 
using a well-defined sample of 81 disk-dominated galaxies from the SDSS redshift survey. We model the 
mass-velocity (TF) and fundamental plane (FP) relations of these galaxies, which are constructed from the 
galaxy stellar mass, disk scale length, and optical Ha rotation velocity at 2.2 scale lengths. We calculate a 
sequence of model galaxy populations, defined by the distribution of the stellar disk-to-total mass fraction, m^. 
We include the effect of adiabatic contraction of dark matter halos in response to condensation of baryons. We 
find that models with constant under-predict the intrinsic scatter of the TF and FP relations and predict an 
(unobserved) strong correlation between TF residuals, even with the full range of halo concentration scatter. 
Introducing a scatter of disk mass fractions and allowing the mean value thd to scale with the stellar surface 
density significantly improves observational match of both the slope and intercept of the model TF relation and 
reduces the predicted residual correlation enough to be statistically consistent with the data. The distribution of 
angular momentum parameters \d required to match the observed disk scale lengths is significantly narrower 
than that predicted for halo spin parameters. Our best-fit models with a Kroupa stellar IMF over-produce the 
galaxy stellar mass function and predict the virial r-band mass-to-light ratios, Mvii/L,-, systematically lower 
than those inferred from galaxy-galaxy weak lensing and satellite dynamics. We investigate three possible 
solutions to these problems: (1) ignoring the effects of adiabatic contraction, (2) adopting a "light" stellar 
IMF with M,/L lower than the Kroupa IMF by 0.15 dex, or (3) considering the lower halo concentrations 
predicted for a low cosmological power spectrum normalization erg w 0.74. In combination with our proposed 
correlation of m,/ with stellar surface density, any of these solutions yields acceptable residual correlations and 
relieves most of the observational tension between the TF relation and the galaxy stellar mass function. 
Subject headings: cosmology: theory — dark matter — galaxies: formation — galaxies: halos 



1. INTRODUCTION 

In the standard theoretical framework of galaxy forma- 
tion, baryons cool, conde nse, and form stars in the cen- 
ters of dark matter halos jWhite &ReeslfT97l . Observed 
disk galaxies obey a tight correlation between luminos- 
ity and rotation s peed, k nown as the Tully-Fisher (here- 
after TF) relation jTullv & Fisher 1977) . The slope, inter- 
cept, and scatter of the TF relation are critical con straints 
on galaxy formation model s (e.g., IC ole & Kaisei* 
Kauffmann et al. 1993; Cole et al.111994 ; Eisenstein & LoeL, 
1996; Somerville & Primack 1 19991 ISteinmetz & Navarro 
1999) . These constraints can be characterized clearly within 
the dissipative co l lapse mod eling framewo rk developed by 
iFall & Efstathiou" Fmd) and "Gunn' ("1983"), updated to the 
cold dark matter (CDM) scenario by Dalcanton et al. ( 1997) 
and Mo et al. ( 1998). In this framework, the disk rotation 
velocity is determined by the ratio, nid, of the disk mass 
to the total halo mass, by the halo density profile, and by 
the angular momentum parameter, A^, which sets the disk 
scale length. In analogy to "fund amental plane" studies of 
the elliptical galaxy population (Diorgovski & Davisj* 119871 
[Dressier et al. 1987), one can consider disk scale length as an 
additi onal parameter in galaxy scaling relations (Shen et al. 
I2002h . Since - GM(R)/R, "maximal disk" models in 
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which baryons dominate the observed rotation curve predict 
a strong anti-correla tion between TF residual and disk size 
(ICourteau & Rix 1999). 

In this paper, we derive empirical constraints on the dis- 
tributions of nid and Xd by modeling the sa mple of disk- 
dominated galaxies from lPizagno et all i2005l) . This sample 
is comprised of 8 1 galaxies with Ha rotation c urves selected 
from the Sloan Digital Sky Survey (SDSS; York & et alj 
EOOO) main galaxy redshift sample ( Strauss & et al. 2002) in 
the absolute magnitude range -18 > Mr > -23, for which 
disk-bulge decomposition yields a best-fit /-band bulge frac- 
tion < 0.1. Pizagno et al. (2005) estimate galaxy stellar 
masses from the observed luminosities and colors, using 
the population synthesis models of Bell et al. (2003) to relate 
the color to the mean stellar mass-to-light ratio. They char- 
acterize rotation velocities by the amplitude of the observed 
rotation curve at 2.2 disk scale lengths, V2.2. We construct 
model galaxy populations with different Piinj) distributions 
and test their ability to reproduce the parameters of the TF 
relation and its residual correlations, after imposing the ob- 
served - Rd distribution as a constraint. Our approach 
differs from previous studies by using a new, homogeneous 
data set with well-defined selection criteria and small obser- 
vational errors, by working with estimated stellar masses in- 
stead of luminosities, and by deriving the P(Xd) distribution 
empirically from the data instead of imposing it a priori from 
theory. While the distribution of the halo angular momentum 
parameter P(X) has been well studied with A^-body simula- 
tions (e.g., Barnes & Efstathiou 1987; Bulloc k et al. 200 laJ, 
the distribution of disk angular momenta could be different 
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because baryons and dark matter exchange angular momen- 
tum, because a biased subset of halo baryons settle into the 
disk, or because disk-dominated galaxies form in a biased 
subset of dark matter halos. An empirical determination of 
PiXd) is therefore a valuable diagnostic of galaxy formation 
physics. 

For our purposes, we define a model of the disk galaxy pop- 
ulation by the probability distribution of m^, which may de- 
pend on stellar mass and scale length, P(md\M^,,Rd). We as- 
sume that the density distribution of dark matter halos is given 
by the NFW profile (Navarro et al. 1997) with a concentration 
parameter c. The distribution function P{c) and mass depen- 
dence of the halo concentrations have been carefully studied 
and are no longer a systematic uncertainty of the models, but 
they influence the scatter and slope of the predicted TP rela- 
tion. For each observed galaxy, we draw random values of 
md and c consistent with the assumed Pinid) and theoretically 
estimated P{c), and calculate the halo mass as M = M^/m^. 
Note that our definition of includes only the stellar mass 
and does not include the contribution of cold gas in the disk. 
We make this choice of necessity because we do not have gas 
mass measurements; in the range of our sample, gas fractions 
are typically ^ 15%. We determine the disk angular momen- 
tum parameter Xd of each galaxy by matching the observed 
disk scale length Rd- 

We judge the acceptability of a model by how well it re- 
produces the observed joint distributions of M»,/?^, and V2 2. 
Since we use the values of stellar masses and disk sizes as 
input to the models, the bivariate M* - Rd distribution is re- 
produced by construction. The slope, amplitude, and scatter 
of the - V2.2 relation and the correlation of residuals of 
the -Rd and M» -V2.2 relations (ARd and AV2.2, respec- 
tively) serve as tests. Alternatively, we consider a "funda- 
mental plane"-type (FP) relation for all three variables, M», 
Rd, and V2.2. After finding the distributions P{md\M^,,Rd) that 
give acceptable fits to the observed relations, we examine the 
ability of these models to fit the observed baryon mass func- 
tion of galaxies. We also make predictions for extended mass 
distributions around the late-type galaxies, which are mea- 
sured by weak lensing and satellite dynamics. 

We include in our models the effect of adiabatic contrac- 
tion (AC) of dark matter halos in response to central conden- 
sation of baryons. The deeper potential well created by the 
cooling and compression of gas attracts more dark matter to 
the central regions of the galaxy and increases its density. In 
the analytic literature on disk galaxy modeling, this effect has 
sometimes been omitted, or its inclusion treated as a "free pa- 
rameter", on the suspicion that the spherical symmetry and 
ad iabatic growth assumptio ns used in the standard derivation 
of iBlumenthal et alJ mM> might make the results invalid in 
the hierarchical merging paradigm of galaxy formation. Re- 
cently, the AC effect has been tested and quantified for galax- 
ies forming by hierarchical merging, using ultrahigh resolu- 
tion gasdynami cs cosmological sim ulations with cooling and 
star formation ( Gnedin et al.'2004*). The effect is not as strong 
as predicted by the original Blumenthal et al. ( 1986) model, 
mainly because the dark matter particle orbits are elongated 
rather circular, but the enhancement of dark matter density in 
the luminous parts of galaxies is robust. The number of galax- 
ies studied at this level of detail remains small, but all avail- 
able analytical and numerical studies show that enhanced den- 
sity of dark matter halos by dissipative baryons is a generic re- 
sult, which is not sensitive to the assumpti on of slow, smooth 
growth of the central galaxy component iTesseit et aljl2()02t 



i Sellwood & M cGau ghll2()()l iChoi et alJl2()()6t iMaccidetal] 
l2006nWeinb erg et al .1 120061) . We are not aware of any hy- 
drodynamic simulations of galaxy formation that do not show 
this effect. 

The inclusion of AC amplifies the effects of disk gravity, 
making it difficult to reproduce the observ ed lac k of corre- 
lation between disk size and rotation speed dCourt eau & Ri^ 
Il999h . By boosting disk rotation speeds relative to halo 
virial velocities, AC also makes it difficult to reconcile 
the observed Tully-Fisher relation with the observed galaxy 
luminosity function, given the halo population of typical 
CDM models; this reconciliation is a long-standing chal- 
lenge to semi-analytic models of galaxy formation (e.g., 
Kauffmann et al. 1993; Cole et al. 1994) because the con- 
flict emerges from a relatively simple, halo-counting argu- 
ment. These observational challenges make it tempting to 
simply omi t AC when creating models of the galaxy popu- 
lation (e.g., ISomerville & Primack|[T999t iDutton et aUlSool 
Given the numerical simulation results cited above, 
however, we regard omitting AC as a radical departure from a 
well established element of galaxy formation physics. It is a 
departure that must be considered if the observations require 
it, but one should consider comparably radical changes (e.g., 
to cosmological parameters or the stellar initial mass function) 
on an equal basis. Throughout this paper, we compare models 
with AC (computed using the modified model of Gnedin et ai] 
120041) and without AC, but we regard the former as "standard 
galaxy formation physics" and the latter as a speculative sce- 
nario. 

2. OBSERVATIONAL RESULTS 

We ana lyze the samp le of 81 disk-dominated galaxies pre- 
sented by P izagno et al. (2005). The galaxie s were selected 
from the main spectroscopic galaxy sample dStrauss & et alJ 
of the SDSS redshift survey and cover a range of ab- 
solute magnitudes -18.5 > Mr > -23.2, with an approxi- 
mately flat distribution. The surface brightness profiles of 
these galaxies in the /- band we re fitted by a sum of an inclined 
exponential disk and a lSersidfl968 ) bulge, and all the galax- 
ies in our sample were chosen to have the disk-to-total lumi- 
nosity ratio > 0.9. This morphological cut, more stringent 
than in most TF studies, allows us to define scale length and 
velocity measures that are insensitive to ambiguities of disk- 
bulge decomposition. It also makes clear the population we 
are modeling: nearly bulgeless galaxies. Rotation curves for 
these galaxies were obtained using long-slit Ha spectroscopy 
on the Calar Alto 3.5-m telescope and the MDM 2.4-m tele- 
scope. 

We take the measurements of stellar mass, disk scale length, 
and rotation speed from Pizagno et al. (2005|. As our mea- 
sure of rotation speed, we adopt the circular velocity at 2.2 
disk scale lengths, V2,2, which according to Courteau ( 199'^ 
produces the tightest TF relation.'* Stellar masses are esti- 
mated from the SDSS /-band Petrosian luminosities and the 
SDSS g-r model colors, both corrected fo r internal extinc- 
tion. We use the models of iBell et al.l d2003h to compute stel- 
lar mass-to-light ratios as a function of g-r color, except that 
we convert from their "diet Salpeter" stellar initial mass func- 
tion (IMF), which is chosen so that stellar disks have t he max- 
imum mass allowed by rotation curve constraints, to a lKroupal 

Our data co ntain small correctio ns to the observational uncertainty of 
V2.2 compared to lPizagno et alJ r200.'?l . as described by Pizagno et al. (2006, 
in prep). 
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Fig. 1 . — Stellar mass - disk scale length relation for the observed sample. 
Dashed lines are displaced vertically from the mean relation (solid line) by 
±1.65 (T, where a is the intrinsic scatter. For a Gaussian distribution around 
the mean relation with standard deviation a, dashed lines would enclose 90% 
of the data points, in the absence of observational errors (which modestly 
increase the dispersion relative to the intrinsic scatter). 

(|2001) IMF, which better represents direct observational esti- 
mates of the IMF. This conversion lowers the stellar mass-to- 
light ratios by 0.15 dex. Thus, throughout most of this paper, 
a galaxy's "stellar mass" is = M»/L,- x L,, where L,- is the 
inclination-corrected Petrosian /-band luminosity and 



logM* /L; = -0.222 + 0.864 (g-r)- 0.15. 



(1) 



We will also consider a model with a "light" IMF, in which 
t he s tellar mass at fixed L,- and g-r is 0.15 dex lower (see 
^3.2> . Galaxy distances, which affect the stellar mass and 
scale length values, are computed from redshifts assuming 
Hq = 70 km s"' Mpc"'. The observational uncertainties in 
M^: and incorporate a Ict distance uncertainty correspond- 
ing to 300 km s"', the typical amplitude of small scale devia- 
tions from Hubble flow (Strauss & Willick 1995). As empha- 
sized by Pizagno et al. (2006) and evident in the plots below, 
an important feature of this sample relative to many previ- 
ous TF samples is that observational uncertainties are small 
compared to the intrinsic scatter, because axis ratios are pre- 
cisely measured and because the minimum redshift cz = 5000 
km s"' ensures that fractional distance uncertainties are small. 
Throughout this paper, we use the notation log for a base- 10 
logarithm and In for a natural logarithm. We use the units of 
solar mass for M,, kpc for Rj, and km s"' for V2.2. 

Figure ^ shows the distribution of sample galaxies in the 
plane of Rd and M,. Figure [2^ shows the V2.2-Af* relation. 
We calculate the parameters of a linear fit 



log V2.2 = a(logM, - logM^o) + b+ N{0, a) , 



(2) 



using the maximum likelihood method described in Ap- 
pendix. The last term in equation (|2j allows for an intrin- 
sic Gaussian scatter of the relationship with standard devi- 
ation a. The offset, logM^o = 10.5, is chosen to minimize 
the error of the intercept, Ab, and to eliminate correlation 
between Aa and Ab. Table [J lists the best-fit parameters: 
a = 0.291 ± 0.013, /7 = 2.23 1 ± 0.006, a = 0.049 ± 0.007. The 



slope and int rinsic scatter a re identical within the errors to 
those given in lPizagno et alJ (12005 ). who used a different im- 
plementation of the same method. Our intercept is different 
because of the slightly different choice for the offset logM»o- 
We have verified that our method reproduces the values b and 
Ab of Pizagno et al. (2005) given their choice of logM^o- 

As Figure[2shows, the disk scale length correlates strongly 
with the stellar mass. The average relation can be determined 
by the same ML method: 



log/?rf = 0.62 + 0.24(logM,-10.5). 



(3) 



We then remove the mean mass dependence of the veloci- 
ties and sizes, and look at the distribution of residuals, AV2.2 
vs. ARd- Figure 0> shows that this distribution is essen- 
tially a scatter plot, with at most a slight positive correlation 
(Spearman's linear correlation coefficient r = 0.23 ±0.14). 
|Q)urteau & Rix ( 1999) emphasized that the lack of correla- 
tion between TF residuals and disk scale length argues against 
the "maximal disk" hypothesis, in which the stellar disk pro- 
vides a large fraction of the rotational support at l.lRj; with 
^ GM^/R, more compact disks would rotate faster at fixed 
M,, creating a strong negative correlation between AV2.2 and 
ARj. The fact that the residual data reveal only a weak posi- 
tive correlation, rather than a substantial negative correlation, 
will play a key role in our conclusions below about the baryon 
fractions of disk galaxies. 

We have also investigated whether a fundamental plane- 
type fit as a function of two variables, and Rj, would lead 
to a smaller intrinsic scatter of ¥2.2- We use our maximum 
likelihood method to fit a linear relation of the form 

logy2.2 = ai(logM,-logM,o)+fl2(log^j-log/;rfo)+^+A^(0,cr). 

(4) 

The error estimates are uncorrected for logM^o = 10.5 and 
log/?,/o = 0.65. The best-fit parameters are given in Table |2] 
ai = 0.271 ± 0.017, 02 = 0.085 ± 0.054, b = 2.233 ± 0.006. 
However, the estimated intrinsic scatter of this relation is es- 
sentially the same as the intrinsic scatter of the TF relation 
above (0.048 vs. 0.049 dex), which supports the notion that 
the TF relation is already a nearly edge-on view of the fun- 
damental plane. The second, radius-related exponent 02 is 
statistically consistent with zero, and it is much smaller than 
the first, mass-related exponent. The relative unimportance of 
the disk size is another manifestation of the lack of residual 
correlation in Figure I2I5, and the sign of the 02 coefficient is 
again opposite to the expectation for self-gravitating disks. As 
a consistency check, if we substitute the mean mass-size rela- 
tion ^ into equation @, we recover the original TF relation 
(eq. 0). 

3. MODELS 

In order to interpret these data, we construct disk galaxy 
models within the framework of a hierarchical ACDM cos- 
mology. We assume that each galaxy resides within a 
dark matter halo whose virial radius, Ryii-, encloses a mean 
overdensity A = 101 relative to the critical matter density, 
appropriate for a flat c osmological model with fim w 0.3 
( Brvan & Normanll998t) . We assume th at the mass profiles o f 
these halos are given by an NFW model dNavarro et alJ 19971) . 
which is parametrized by the virial mass, M, and concentra- 
tion parameter, c. Halos found in dissipationless ACDM sim- 
ulations have concentrations whose mean c correlates with 
halo mass and whose distribution has a Gaussian scatter in 
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Fig. 2. — Left panel: Observed TF relation in the i band with the lBell et all 12003) mass-to-Ught ratios modified to a Kroupa IMF. Dashed lines are displaced 
vertically from the mean relation by ±1.65(j, where a is the intrinsic scatter. For a Gaussian distribution with standai'd deviation cr, the two dashed lines would 
enclose 90% of data points, in the absence of observational errors. Right panel: Residuals of the TF relation, A V2.2, vs. residuals of the mass-size relation, AS^. 
The dashed ellipse is oriented in the direction of the residual correlation, with its major and minor axes equal to 1 .65 times the rms deviation in the parallel and 
perpendicular directions, respectively. This ellipse would enclose 81% of data points if the distribution was Gaussian. 



logc: 

P{c)d\ogc= exp — d\ogc. (5) 

VIttUc \ 20-2 J 

iBullock etaH (1200 Ibh find an empirical fit c(M) w 
12.8(M/10'2r'MQ)-*'-i3 for N-hody simulations with 
rim = 1 - f^A = 0.3 and ag = 1 .0. We adopt this c{M) relation 
as our standard, but in ^we also consider the lower con- 
centrations predicted by the Bullock et al. (2001b) analytical 
model in the case of lower a%. The scatter is CTc = 0.14 dex 
for all halos, and may be smalle r for those halos tha t have 
not had a recent major merger (IWechs ler et al."2002). As 
discussed in the introduction, our standard models include 
the effect of contraction of dark matter hal os in response to 
condensation of baryons in their centers ([ Barnes & White! 
[1^84; Blumenthal et al. 1986; Rv den & Ciunnlll987h . cam- 
puted using the modified adiabatic contraction model of 
[Gnedin etal. (2004). We also explore models in which 
adiabatic contraction does not occur, although we consider 
these less physically reasonable. 

We treat the stellar component of each galaxy as an expo- 
nential disk with the measured mass and scale length R^; 
since disk-bulge decomposition yields a bulge-to-total ratio 
< 0.1 for all galaxies in our sample, this assumption is appro- 
priate. We do not have gas masses available for our sample, 
and therefore we omit the gas contribution to the disk poten- 
tial and denote the ratio of the disk stellar mass to the halo 
virial mass by nid- In many theoretical papers, nid represents 
the fraction of disk baryon mass (stars plus cold gas). Since 
gas fractions are typically small in the luminosity range that 
we consider, the difference between th e two definitions sh ould 
only be 15% on average (based on the lKannappanl2004l rela- 
tion for u-r galaxy colors). 

The choice of and c completely specifies the model of 
an individual observed galaxy, given M, and The angular 
momentum parameter \d is not an input, but can instead be 



derived from the total disk angular momentum: 

Jd = \d VlM^RvuVv,, = J VAR)R2TTT.(R)RdR, (6) 

where S(/?) is the surface mass density of the disk at az- 
imuthal radius R, and Vc(R) is the rotation velocity in the plane 
of the disk, which includes the contributions of the flat disk 
and the spherical halo. Note that we use the model Vc(R) in 
computing from equation (|5}, and a surface density profile 
T^iR) that is exponential with the observed M^, and Rj. The 
virial velocity is defined as V^^^. = GM/R^i^, where /?vir is the 
virial radius. The angular momentum of the halo can be sim- 
ilarly expressed using the halo spin parameter A: 

Ji, = XV2MRvnV,u. (7) 
In our framework, a model of the disk galaxy population 
is defined by the probability distribution P(m,/|M*,7?j), and 
our goal is to investigate what forms of this distribution yield 
an acceptable fit to the observed distribution of galaxies in 
the (M*,/?(/, V2.2) space. For each individual galaxy of the 
IPizagno et a l. (2005) sample, we proceed as follows. (1) 
We draw a disk fraction, nid, from the assumed distribution 
P{md\M^^Rd) and calculate the halo virial mass, M = M^,/md, 
using the observed stellar mass M*. (2) We draw an NFW 
concentration parameter, c, from the theoretical distribution 
Pic), given by equation (|5}. In some models we set the scatter 
of concentrations to cTc = 0. (3) We calculate the adiabatically- 
contracted halo mass profile using the modified model of AC 
with the known M*, M, c, and Rd- We also run models omit- 
ting the AC effect. (4) We calculate the rotation velocity at 
2.2Rd and add random measurement errors to logM,^, \ogRd, 
and logV2.2, drawn from their respective Gaussian distribu- 
tions, each with the standard deviation equal to the corre- 
sponding uncertainty for the observed galaxy. This step is 
necessary in order to correctly model the correlation of resid- 
uals. 

In constructing our model sample, we implicitly assume 
that the model predictions can be calculated with sufficient 
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Fig. 3. — TF relation and residuals for three representative models with constant mj and no scatter in halo concentrations, ac = 0. Filled circles ai'e for models 
including the effect of adiabatic contraction of dark matter, open squares are for models without the AC effect. Dashed lines (left panel) and contours (right 
panel) are for the data from Fig. |2| Scatter in the light panels is almost entirely caused by the measurement en'ors added to model values of M», Rj, and ¥2,2- 
The solid lines in the right panels show the original model realizations without including the measurement errors. These lines coiTespond to the case with AC, 
the case without AC is similar 



accuracy and that the only uncertainty in comparing with the 
data is due to the measurement errors. However, some of our 
models are probabilistic in nature, as they contain scatter of 
halo concentrations and disk fractions. Given only 8 1 galaxies 
in the observed sample, the predictions of these models vary 
from one realization to another. Even for fixed nid and c(M), 
the random errors added to logM* and logV2.2 lead to small 
variations in best-fit parameters among realizations. In order 
to suppress such statistical variations and obtain robust aver- 
age model predictions, we take a large number of repetitions 
of each observed galaxy, A'rep, and fit a single TF relation (and 
FP relation) to the combined sample of 8 1 x A',ep galaxies. We 
have conducted a convergence study, by increasing A'lep until 
the model fit parameters vary by less than 10% of the corre- 
sponding observational uncertainty. We have found that this 
condition requires at least 100 realizations of the sample in 
some cases, and we adopt A^,ep = 100 for all our models. 

We fit the TF relation to our sample of 8 100 model galaxies 
as in equation (|2j, with the same logM^o used for the obser- 
vational data. We compute AV2.2 and ARj residuals relative 
to this best-fit TF relation and compute the correlation coeffi- 



cient as for the observational data. The best-fit slopes a, zero- 
points b, intrinsic dispersions a, and correlation coefficients r 
of all our model samples are given in Tabled We also fit the 
FP relation similarly and list the best-fit parameters in Table 

m 

For a given model sample, we evaluate the goodness-of-fit 
to the observed TF relation based on the statistic: 

2 / ^mod ^obs \ / ^mod ^obs \ / ^mod <^obs 



X = 



\ Aflobs 
/ ^mod ^obs 



(8) 



V Arobs 

For the goodness-of-fit of the FP relation, we use 

2 / ^l.mod ^l.obs \ , / ^2, mod ^2.obs \ , / ^mod ^obs 



X = 



Afli.obs 



f^mod *^obs 



Afl2,obs 



V A/7obs 



V AfJobs 

Note that all of our models reproduce the observed (M* , Rj) 
distribution by construction. 
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Fig. 4. — Best-fit models to the TF relation and residuals with constant mj, including the scatter of halo concentrations, (Jc = 0.14. Bottom panel: models 
including the effect of halo contraction, with nv = 0.08; top panels: without the AC effect, with mj = 0.04. 



3.1. Constant disk mass fraction 

First, we explore models with a constant disk fraction 
and no scatter of halo concentrations, <Jc = 0. In this case the 
model predictions for each galaxy are deterministic, except 
for the added measurement errors. Figure |3] shows the TF 
relation for model galaxies and the distribution of residuals. 
The top panels correspond to a low disk mass fraction, = 
0.02. For a given galaxy stellar mass, low nid corresponds to a 
large dark matter halo and, therefore, a high circular velocity. 
The model velocities lie consistently above the observed mean 
relation, indicating that the true average value of nid is larger 
than 0.02. Since the halo dominates the mass even at small 
radii, the adiabatic contraction effect is not very important, 
and the model realizations with and without the AC effect are 
similar to each other. 

Since the disk contribution to the rotation curve is small, the 
disk size has little impact on the velocity residuals (top right 
panel of Fig.|3}- Most of the plotted velocity residuals arise 
from the random "measurement" errors added to the model 
galaxy values. A solid line in that panel shows the original 
predicted relation of the residuals, before we add the random 
errors. If the disk completely dominated the rotation curve, at 
a given stellar mass we would expect an anti-correlation of the 
form AlogV2.2 = -0.5 Alog/?^/. With a large contribution of 
the dark halo, the relation of the residuals is much flatter but 
is still approximately monotonic. The relation plotted by the 
solid line is insensitive to the stellar mass - we have split the 
sample in bins of large and small and found approximately 
the same relation as for the whole sample. 

In statistical terms, the = 0.02 models have accept- 
able TF slopes but zero-points that are too high, b = 2.306 
(2.265) with (without) AC, compared to the observed values 
of 2.231 ± 0.006. Furthermore, both models have an intrinsic 
scatter that is too small, a = 0.014 (0.009) compared to the 



observed a = 0.049 ± 0.007. 

The middle and bottom rows of Fig.|3]show model realiza- 
tions for larger disk mass fractions, = 0.04 and = 0.08, 
respectively. Increasing lowers the halo mass, given the 
observed M», and thus lowers the TF zero-point, bringing the 
models into better accord with the data. 

Within this class of constant-m^/, constant-c models, the 
best fit is achieved for ~ 0.08 with AC and k, 0.04 
without; AC boosts V2.2 and therefore requires higher m^. 
However, both models are statistically unacceptable, with to- 
tal of 36 and 46, respectively (for four observables and 
one adjustable parameter, hence three degrees of freedom). 
The no-AC model has a slope that is slightly too steep (0.295 
vs. 0.291 ±0.013), but its more serious problems are pre- 
dicting too Uttle scatter (a = 0.016 vs. 0.049 ±0.007) and 
predicting anti-correlated residuals (coefficient r = -0.26 vs. 
+0.23 ± 0.14). Disk gravity plays a much larger role in the 
AC model, because of both the higher and the effects of 
AC itself. The slope discrepancy becomes worse (a = 0.322) 
because more massive disks are higher surface density on av- 
erage and produce stronger boosts in ¥2.2- The intrinsic scat- 
ter is closer to the observed value (cr = 0.034), but this scat- 
ter is driven entirely by variations in disk size, so the anti- 
correlation of residuals is worse (r = -0.57). 

The conclusions of the FP analysis are similar. The mass 
coefficient oi increases and the intercept b decreases with in- 
creasing nid- The radius coefficient a2 becomes more nega- 
tive, further away from the observed value, the signature of 
an overly dominant disk. The intrinsic scatter of the FP re- 
lation does not vary noticeably with mj, but it is much lower 
than observed, a « 0.008 vs. 0.048 ± 0.006 (see TableEJ. 

It is easy to see why models with fixed values of and 
c cannot reproduce the observed intrinsic scatter and uncor- 
related residuals, regardless of the value of and the pres- 
ence or absence of AC. In these models, the pre-contraction 



TuUy-Fisher Relation 



7 



halo profile is determined by M*, and the variation in is 
therefore the only available source of intrinsic scatter of V2.2 
at fixed M*. Thus, if there is significant intrinsic scatter, there 
must be strong TF residual anti-correlations. The FP analysis 
removes the scatter associated with R4, making a independent 
of ifid (and much too small). 

Of course, dark matter halos in cosmological simulations 
display a range of concentrations at fixed virial mass, char- 
acterized by the log-normal distribution of equation (|5} with 
(Tc = 0. 14 dex. We add this variation to our second set of mod- 
els, plotted in Figure 0] The scatter of halo profiles signif- 
icantly increases the intrinsic scatter of the model TF rela- 
tion (see Table[0, especially for low m,/ models, and brings it 
within striking range of the observed value. Since the concen- 
tration variations are uncorrected with Rd, they also reduce 
the anti-correlation of TF residuals. Both effects improve x^, 
but they still do not produce acceptable fits to the data, for ei- 
ther the TF or FP relations. The best models have total x' of 
40 with AC and 13 without (Tabled TF) and larger values for 
FP (TableEJ, again for three degrees of freedom. 

Note that in the Mo et al. (1998) formalism, the concen- 
tration parameter encodes the effect of halo formation his- 
tory on halo profile, essentially the ef fect d escribed as for- 
mation redshifi in earlie r work (e.g.. ICole & Kaiser 1989; 
lEisenstein & Loe? While these earlier papers high- 

lighted the problem of over-predicting the scatter of the TF 
relation, because of variations in formation redshift, the range 
of concentrations found in cosmological A^-body simulations 
is actually insufficient in itself to explain the observed scatter 
Furthermore, since the tail of low concentrations comes from 
halos with active late-time merger histories, we consider it 
plausible that the concentration parameters of halos that host 
disk-dominated galaxies have smaller scatter than CTc = 0.14, 
perhaps much smaller 

Suppose that we have systematically underestimated the 
"observational" errors in V22 and M*, perhaps because of 
non-circula r motions, disk ellipticitie s that affect inclination 
corrections jFranx & de Zeeuw|[T99l . or variations in stellar 
m ass-to-light ratio at fixed color (present at the 0. 1 dex level 
in lBell et al.l2003l) . We put "observational" in quotes because 
these could equally well be considered "physical" sources of 
scatter, but they are in any event not included in the models 
above. We consider an upper limit on these effects by adding 
random velocity errors to the models until the intrinsic scatter 
of the model relationship matches that of the observed TF re- 
lation. In practice, we add log-normally distributed errors in 
V2.2 with the variance equal to al^^ - cr^^^- The TF intrinsic 
scatter matches the observed value by construction, and the 
FP scatter is close to that observed. However, the TF residual 
anti-correlations remain much too strong, and the ai slopes of 
the FP fits remain negative and large (Tables ^ and |2}- Re- 
alizations of these models (not shown) are visually similar to 
Figure 0] 

We conclude that models with a single value of cannot 
reproduce the observed TF relation and TF residuals of our 
disk-dominated galaxy sample, and that much of the intrinsic 
scatter of the TF relation must arise from scatter in the ratio of 
stellar disk mass to halo mass. Pizagno et al. ( 2005) reached 
the same conclusion based on a more qualitative comparison 
of models to the data. We turn to models with m^z-scatter 
in the next section, but first we examine the distributions 
inferred for our constant m,/ models. 

Figure |5] shows the P{\d) distributions for the two models 
illustrated in Figure^] which include concentration scatter In 
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Fig . 5 . — Distribution of disk angular momentum parameter required by 
our models to reproduce the observed disk sizes for the models with constant 
m^, including scatter of halo concentrations, ac = 0.14. Solid histogram is 
for nv = 0-08 with AC, dashed histogram is for mj = 0.04 without AC. For 
comparison, dotted hatched histogram shows the P(X) distribution of dark 
matter halos in dissipationles s co smological simulations, parametrized by a 
log-normal function (see eq. 1101 in the text). Normalization is such that the 
peak of the first distribution is unity. 

the nid = 0.08, AC model, galaxies reside in less massive ha- 
los with smaller virial radii, so larger values are required 
to produce their observed scale lengths. If the models ex- 
actly reproduced each galaxy's observed rotation curve VdR), 

then equation (|6|l implies oc R~l^V~:^ oc nij^ on a galaxy- 
by-galaxy basis, independent of AC except for its impact on 
the best-fitting nid- Since the two models approximately re- 
produce the observed TF zero-point and slope, they obey this 
scaling approximately but not exactly. The modes of the in- 
ferred Xd distributions are = 0.056 (AC, nid = 0.08) and 
Xd = 0.036 (no AC, nid = 0.04), and the standard deviations of 
\nXd are a\d = 0.30 and a\d = 0.36, respectively. (In agree- 
ment with established notation, we quote here the standard 
deviation of the natural logarithm of A; all the other disper- 
sions in this paper are for base-10 logarithms.) 

The shaded histogram in Figure |5] shows the theoretically 
predicted distribution of dark matter halo spin parameters (eq. 
171 ) derived from cosmological A^-body simulations (e.g., 
.Vitvit ska et al. 2002), 

, , 1 / [lnA-lnA]2\ c/A 

PiX)dX=^^cxp{- ' M — , (10) 

with A = 0.0 45, ax = 0.56. According t o tidal tor que theory 
(e.g., l Hovlelf!94l lPeebles|[T969: Efstathiou & Jones 198Q) . 
baryons should acquire the same specific angular momentum 
as dark matter at the time of halo turnaround (when most of 
the angular momentum is imparted). However, as noted in 
the introduction, the Xd distribution can depart from the A 
distribution because disks acquire a biased subset of the halo 
baryons, because disk galaxies form in a biased subset of dark 
matter halos, because angular momentum is transferred be- 
tween the dark matter and the baryons, or perhaps because 
assembly along a filament brings baryons from beyond the 
virial radius. Figure |5l suggests that one or more of these ef- 
fects must play an important role in determining disk angular 
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Fig. 6. — Best-fit models to the TF relation and residuals with the scatter of disk mass fractions, cr,,, = 0.22 (top panel, no AC) and cr,,, = 0.27 {bottom panel, 
with AC), and scatter of halo concentrations, = 0.14. 



momenta. In both of our models, the inferred distribu- 
tions are narrower than the theoretical A distribution. Interest- 
ingly, Wec hsler et al. (2002) and D'Onghia & Burkert (2004) 
find in their cosmological simulations that halos that have not 
experienced a recent major merger, and therefore are likely 
to host disk-dominated galaxies, have a narrower distribution 
P(X) than all halos. The absence of a low-A^ tail can plausibly 
be explained by low-spin halos producing galaxies with sub- 
stantial bulges. For the AC model with = 0.08, the inferred 
Xd values are systematically higher than the A^-body A values. 
Since angular momentum transfer seems more likely to re- 
duce Xd below A (Navarro & Steinmetz 1997), this difference 
is more plausibly explained by the other effects mentioned 
above. 

3.2. Log-normal distributions of disk mass fractions 

In the lMo et alJ ((193 8|) formalism of disk galaxy modeling, 
the sources of scatter in the TF relation are halo concentra- 
tion, disk size, and disk mass fraction. We have shown in 
the previous section that the first two of these predict too lit- 
tle scatter and too much residual correlation, compared with 
the data. Even though boosting the "measurement" errors of 
model velocities, beyond the level expected in the data, can in 
principle explain the observed scatter, it does not remove the 
residual correlation. In this section, we consider models with 
a log-normal scatter cr,,,. 



1 



27rcr„ 



[logm^-logm^] 



P{md)dlogmd = —f= — exp ( —^^ j dlogmd, 

(11) 

and allow nid to scale with stellar mass or stellar surface den- 
sity. In all cases with nid scatter, we truncate Pimd) above the 
adopted maximum value of 0. 15 (close to the universal baryon 
fraction) and below the minimum value of 0.001. 
We first fix rhd at the previous best-fit values for the cases 



with and without AC, and increase the dispersion until we 
match the observed intrinsic scatter of the TF relation. This 
match requires a,,, w 0.2-0.3. Figure |S] shows that even in 
this case the residual anti-correlation persists, especially in the 
realization with AC. Similarly, in the FP analysis the slope 02 
remains consistently negative, instead of positive as required 
by the data. 

In addition, the AC model has a significantly discrepant 
TF slope (fl = 0.320 instead of 0.291 ±0.013) and intercept 
(Z7 = 2.249 instead of 2.231 ±0.006). These are 2-3 ct statis- 
tical deviations. The steeper slope is visually evident in the 
bottom left panel of Figure |6l The model and observed re- 
lations cross, but at a mass below M^o where the intercept is 
defined. Relative to the constant-m^ models, the total de- 
clines from 40 to 35 in the AC model and from 13 to 10 in the 
no- AC model, but since adding a„, as an adjustable parame- 
ter decreases the number of degrees of freedom from three to 
two, the reduced actually goes up in each case. 

The efficiency of gas cooling in galactic halos or gas ejec- 
tion or heating by feedback could plausibly depend on galaxy 
mass and density. Accordingly, the typical disk mass fraction 
may vary with these parameters, and Pizaano et al. (2005) ar- 
gued for lower average m^ values at low stellar masses. Next, 
we consider a model where the central value rhd scales with 
the stellar mass as 



md = nid.o 



10'"5M, 







(12) 



We have run a grid of models, varying the normalization ifid.o, 
the log-normal width cr,,,, and the exponent pi, to find the 
minimum x^- Figure shows the best-fit models for the 
cases with and without the AC effect, for which the exponents 
Pi = 0.3 and pi = 0.1, respectively. By boosting the circular 
velocities of low M^, galaxies, the scaling with stellar mass 
significantly improves the slope and normalization of the TF 



TuUy-Fisher Relation 



9 



500 



0} 
> 



100 

50 
500 



TTTTj 1 — I I I I lllj 1 I I I I lllj 1 I l>llll 

m, = 0,04 .-^ 
CT„=0.14 □ ^iftrfi/ 



yd 

id I I 



I 



100 
50 



TTTTl 1 — I I I I III 



m, = 0.1 
a,=0.14 
AC 



I I 1 1 iiiij — I I I >i 1 1 1 



LlJ I I 1^ I ml 



I 



> 



> 
o 



109 lo'o 10" 
M. (MJ 



10 



12 



0.2 
0.1 



-0.1 

-0.2 
0.2 

0.1 



-0.1 
-0.2 



I ' ' ' I 



I I I I I I I 



I ' ' ' I 



' [IE 



° Di, □ / 



I I I I I 



I I I I I I I 



I I I I I 



I ' ' ' I 



I I I I I I I 



I ' ' ' I 



4 . •* *• t . . ^ 



1 1 1 1 1 1 1 1 1 1 1 



-0.4-0.2 0.2 0.4 
AlogR, 



Fig. 7. — Best-fit models to the TF relation and residuals with the mean disk mass fraction increasing with stellar mass, oc Mj' , including the scatter a,,, 
and (Tc = 0.14. Top panel: pi = 0.1, a„, = 0.25, no AC; bottom panel: pi = 0.3, a,,, = 0.21, with AC. 



relation for models with AC (both are now within 1 cr of the 
data). However, the predicted anti-correlation of the residuals 
is still inconsistent with the data, as is evident from the lower 
right panel of Figure^ The no-AC model with M» scaling is 
consistent with the observed slope, intercept, and scatter, but 
it is still marginally inconsistent with the residual correlation 
(r = -0.17 vs. 0.23 ±0.14). 

These models now have three adjustable parameters {mdfl, 
pi, am, with tJc fixed by theory) and the degeneracy among 
the best-fit values is substantial. Within the limit Ax^ = 1, 
the range of allowed model parameters is Ap; « 0.1, Act,,, w 
0.2-0.3. Our inferred parameters should therefore be taken 
as illustrating a preferred, but by no means unique, trend. 
Note that even three adjustable parameters do not allow us 
to achieve an acceptable match to the four TF parameters (a, 
b, CT, r) or to the four FP parameters. Roughly speaking, m^.o 
affects the intercept (high nid - low b) and residual correlation 
(high nid - strongly anti-correlated residuals), ct,„ affects the 
scatter (high ct^ - high scatter), and pi affects the slope (high 
pi - low a), with CTc producing additional scatter, but no ele- 
ment in these models can remove the residual anti-correlation. 

Disk galaxy models of Firmani & Avila-Reese (2000) show 
a trend of increasing disk contribution to the rotation curve 
with increasing stellar surface density, which in their mod- 
els arises from the correlation of the star formati on efficiency 
with the surface density. In observational studies. lZavala et aP 
(|2003) and Pizagno et al. (2005) find that the ratio of halo 
mass to stellar mass within the luminous regions of disks cor- 
relates strongly with disk surface density. Motivated in part 
by these results, we consider a set of models where the central 
value ifid scales with the stellar surface density, E» = M^.R~^: 

_ _ ( M^Rf 

V 10^-2 Mokpc ^/ 
The best-fit scalings for p2 are 0.65 (AC) and 0.2 (no AC), 



respectively. Figure|8]shows that such models reproduce both 
the TF and FP relations fairly well, even including the AC ef- 
fect. Given the mass-size relation (eq. 0), distributions il2\ 

and (I13t are internally consistent; oc M^^ and p2 ~ 2pi 
(see Table [Q. However, tying irid to E, leads to substan- 
tially better agreement with the data, greatly reducing the 
anti-correlation of TF residuals. The model without AC is 
consistent with the observed slope, intercept, and scatter at 
the 1 CT level, while its predicted residual anti-correlation re- 
mains ~ 2ct discrepant, r = -0.07 vs. 0.23 ± 0.14. The FP 
parameters are consistent at the 1 ct level, except for the ra- 
dius coefficient, which is a2 = 0.022 vs. 0.085 ± 0.054. The 
AC model remains notably discrepant with the TF residual 
correlation, r = -0.18, and the FP radius scaling, 02 = -0.022, 
but these discrepancies are much smaller than in any of our 
previous models, and other parameters agree within 1 ct. If 
we minimize with respect to the FP parameters instead of 
the TF parameters, we obtain somewhat stronger density scal- 
ings; these higher p2 models are also listed in Tables^andl^] 
Surface density scaling improves the TF and FP fits because 
it places higher surface density galaxies in lower mass ha- 
los. The anti-correlation between disk-induced rotation ve- 
locity and halo-induced rotation velocity counters the trend 
for more compact disks to rotate faster - in effect, there is 
no correlation between radius residual and velocity residual 
because two correlations of opposite sign cancel each other 
While this balance of competing effects may seem coinciden- 
tal, a correlation of ifid with E* seems eminently reasonable 
on physical grounds. Cooling and star formation should be 
more efficient for denser baryon distributions, and supernova 
feedback is probably less efficient at ejecting gas from higher 
surface density disks. In fact, the causal arrow may point in 
the other direction - high E,^ disks may have higher surface 
density largely because they have higher nid- Finally, since 
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Fig. 8. — Best-fit models to the TF relation and residuals with the mean disk mass fraction increasing with stellar surface density, mj oc T,p, including the 
scatter a„ and cTc = 0.14. Top panel: p2 = 0.2, it„, = 0.25, no AC; bottom panel: P2 = 0.65, cr,„ = 0.2, with AC. 
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Fig. 9. — Best-fit models to the TF relation and residuals for the case of light IMF, with the mean disk mass fraction increasing with stellar surface density, 
rSrf oc S^**, including the scatter cr,„ = 0.25 and ac = 0.14, and AC. The dashed lines representing the observations in the left panel ai'e shifted by -0.15 dex in 
M« to account for the changed IMF. 



we define in terms of disk stellar mass, a trend of m^/ will 
arise automatically from the larger gas fractions of low sur- 
face brightness galaxies. This trend is driven physically by 
the correlation of star formati on rate with gas surface density 
(ISchmidt 1959; Kennicutt 1998). 

Models without AC constitute a variation of our standard 
models (with AC). As another variation, we consider models 
with a light stellar IMF. We reduce stellar mass-to-light ratios 
(eq. 10) by 0.15 dex relative to the Kroupa values. The less 
massive stellar disks now require more massive dark matter 
halos to reproduce the observed V2.2, and because adiabatic 
contraction has less impact, the best-fit values of nid decrease 
by a factor substantially larger than 10"^^. Figure |9] shows 
the best-fit model with AC and w,/ cx E^ ^^. The normaliza- 
tion in equation il3\ drops from fhjx) = 0.1 to mjx) = 0.04. 



This model has no TF residual anti-correlation and the right 
amount of intrinsic scatter, although the slope and intercept of 
the TF relation depart from the observed values more than in 
our standard Kroupa-IMF model. This model is also statisti- 
cally consistent with the observed FP relation, reproducing its 
mass coefficient, intercept, and scatter, and lying within 1. 2 a 
of the radius coefficient 02- Although we have no additional 
evidence to justify using the reduced mass-to-light ratios, this 
alternative provides a good description of our galaxy sample. 
If the population synthesis models themselves are slightly in- 
accurate, in the sense of predicting M»/L,- values that are too 
high for a given g-r color, then part of the desired 0.15-dex 
reduction might arise even without a change to the IMF. 

One of the persistent issues in disk galaxy studies is the role 
of disk gravity in determining the rotation curve. In "maximal 
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disk" models, the disk typicall y contributes 75% -95% of the 
rotational support at 2.2Rd ( Sackettl[l997h . Solid points in 
the left panel of Figure [TO] show the ratio of the circular ve- 
locity from the stellar disk to the measured circular velocity 
V2.2 for our obser ved sample, assumin g a Kroupa IMF (sim- 
ilar to Fig. 7 of iPizagno et alJl2005l) . The contribution of 
stars to the circular velocity at 2.27?^ increases monotonically 
with the stellar surface density. Most of the galaxies with 
> 3 X lO^M0kpc~^ are maximal disk by Sackett's crite- 
rion, while the less dense disks are dark matter-dominated.^ 
Open points show the same ratio for the model realization 
shown in the lower panel of Figure |8] with rhd cx S'''^^ and 
AC. The model shows the same trend as the data, except that 
there are essentially no model galaxies with Vs/V2.2 > 0.8. Ra- 
tios higher than this imply such small halos (relative to the 
galaxy baryon mass) that our model never predicts them. The 
right panel of Figure^] shows the same results for the light 
IMF. Note that the observational data points also change in 
this case, because the value of inferred for each galaxy is 
lower and the values of and are therefore lower as well. 
Only two observed galaxies have Ysl^ii > 0.8 in this case, 
while the rest have sub-maximal disks. The correlation of the 
ratio Vs/V2.2 with is stronger than with or other avail- 
able galaxy parameters. This reinforces the assumption of our 
models that the surface density is a major factor regulating the 
stellar mass fraction of disk galaxies. 

Figure shows the distributions of the disk angular mo- 
mentum parameter for the two models of Figure |8l These 
histograms are qualitatively similar to those for the constant- 
nid models shown in Figure |5] The typical values again 
increase with the typical values of m^. _ The mode values 
are Ad ~ 0.054 for the AC model and Xd w 0.036 for the 
no-AC model, compared to A 0.045 for dark matter ha- 
los. For visual clarity, we have omitted the light-IMF model 
from Figure ^2 but its A^ histogram is intermediate between 

^ When using the number 3 X lO' Mq kpc"^, one should keep in mind our 
definition of the surface density, E, = M* R^^, not containing any factors of 
2n. 



the other two, with Aj ~ 0.040. All three model distribu- 
tions are narrower than the corresponding P{X) for dark mat- 
ter halos: axd ~ 0.25 (AC), 0.39 (no AC), 0.34 (light IMF) 
vs. a\ w 0.56. Furthermore, our procedure should overesti- 
mate the true width of the A^ histogram because we choose a 
random value of nid from P(md\M^,,Rd), thus adding artificial 
scatter to the Xd values given the Rd constraints. (We do not 
expect to draw the correct value of nid on a galaxy-by-galaxy 
basis, even if our model incorporates the correct spread of nid 
values.) To test the magnitude of this effect, we have applied 
our procedure to a realization of the AC model, choosing new, 
random values of irid and c for each model galaxy. The in- 
put Xd distribution has a standard deviation a\d = 0.25, while 
the distribution inferred by fitting the model data points has 
<T\d = 0.34, nearly 40% larger. The mean value A^ is repro- 
duced correctly. Thus the A^ distributions shown in FigureFTTI 
place strong constraints on the allowed values of the angular 
momentum parameter of disk-dominated galaxies. These dis- 
tributions are quite different from the A distribution of A^-body 
halos, implying that only a subset of halos form disk galaxies 
like those in our sample, or that disk baryons have systemati- 
cally different angular momenta from the dark matter. 

Button etal. (2006, hereafter DVDC) have recently car- 
ried out an investigation that is similar to ours in many re- 
spects, though in spirit theirs is more a theoretical modeling 
study constrained by data while ours is an empirical inference 
study informed by theory. The observational results from their 
large compilation of three /-band TF data sets (~ 1300 galax- 
ies with Ha rotation curves) agree well with those from our 
smaller, homogeneous data set, which is reassuring given the 
numerous differences in sample selection, observations, anal- 
ysis methods, and typical galaxy distances (which strongly af- 
fect the observational error budget because of peculiar veloc- 
ities). In particular (see Fig. 3 of DVDC) we find similar TF 
slope and intercept and similar lack of correlation between TF 
residual and radius residual; our galaxies have slightly larger 
average scale lengths at fixed luminosity, which is plausibly 
an effect of our stricter limit on bulge-to-disk ratio. DVDC's 
estimated intrinsic scatter of 0.052 dex in logV2.2 at fixed /- 
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Fig. 1 1 . — Distribution of disk angular momentum parameter required by 
our best-fit models with oc E','" with AC (solid histogram) and without AC 
(dashed histogram). Dotted hatched histogram shows the P(X) distribution of 
dark matter halos as in Fig. |5] 



band luminosity, inferred by subtracting observational errors 
in quadrature from the total scatter, is only marginally larger 
than our maximum-likelihood estimate of the intrinsic scatter 
at fixed stellar mass, a = 0.049 ± 0.007 dex. This similar- 
ity is somewhat surprising, since the scatter at fixed luminos- 
ity should be increased by galaxy-to-galaxy scatter in /-band 
mass-to-light ratios. Furthermore, if DVDC included distance 
uncertainties in their observational error budget, their estimate 
of the intrinsic scatter would decrease, coming closer to our 
value or perhaps even falling below it. 

There are many differences of detail in our theoretical mod- 
eling methods: we model relations involving estimated stellar 
masses while DVDC use /-band luminosities; we infer dis- 
tributions by fitting individual values while DVDC assume 
a log-normal distribution of adjustable width; and we adopt 
a lKrouDal ( l200lt) IMF as our standard while DVDC take Bell 
& de Jong's (2001) "diet Salpeter" IMF as fiducial. Despite 
these differences, our qualitative theoretical conclusions are 
largely concordant to the extent that they can be compared. 
Both studies infer a narrow distribution of relative to N- 
body predictions for the halo A distribution. DVDC do not 
consider models with nid connected to S*. However, their 
surface density threshold for star formation has an effect of 
the same sign (since they define in terms of baryonic mass 
while we define it in terms of stellar mass), and they find that 
this threshold plays an important role in suppressing TF resid- 
ual correlations. We find that scatter in nid is required to re- 
produce the observed TF scatter, while DVDC do not. Given 
the similar observational estimates of the intrinsic scatter, this 
difference probably arises largely because DVDC incorporate 
scatter in /-band stellar mass-to-light ratios while we use g-r 
color to estimate M/L, and assume no scatter about the mean 
relation. DVDC emphasize the galaxy luminosity function 
as an additional constraint on models, and we will take up 
this point ourselves in l^below (see Figure [141 . We concur 

^ As we were already finalizing our manuscript when the DVDC preprint 
became available, we have not carried out a comprehensive comparison of 
results. 



with DVDC that omitting adiabatic contraction, adopting a 
lighter IMF, or lowering halo concentrations (see all help 
to improve agreement with the observed TF relation (includ- 
ing residual correlations) and to reduce the tension between 
the TF relation and the galaxy stellar mass function (or lu- 
minosity function). We are more skeptical about the physical 
plausibility of omitting (or, in DVDC's favored model, revers- 
ing) the effects of adiabatic contraction, for the reasons cited 
in ^ so we lean towards some combination of light IMF and 
lower halo concentration as a more likely explanation. 

3.3. Individual galaxy ( "exact") fit models 

A common approach in modeUng the TF relation has been 
to choose a distribution P( Xd) a pri ori ( Kauffmann et al 
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iDalcanton et al .''1997; Mo et al. 1998; vandenBoscl 
In the previous sections, we have modified this ap- 
proach by choosing Xd to match the observed disk scale length 
Rd on a galaxy-by-galaxy basis. In doing so, we have derived 
an empirical estimate of P(Xd), given the model distribution 
P(md\M^,Rd). We can go one step further and choose for 
each galaxy to exactly reproduce its observed velocity ¥2.2- 
We still assume that the dark matter halos are described by 
the NFW model with a concentration parameter dependent on 
M. Models constructed in this way reproduce the observed 
M^,, Rd, and V2.2 distributions by construction, and they yield 
P{md) as a result. 

Jimenez et al. (2003) follow this approach to estimate 
the disk mass fractions and angular momentum parame- 
ters from a large sample of Ha rotation curves, assembled 
from various sources. They consider cuspy NFW mod- 
els and cored pseudo-isothermal spheres for the dark halos 
and apply the st andard model of adiabatic contraction (e.g., 
iBlumenthal et alJil98 6). In contrast to our procedure, they do 
not use stellar mass-to-light ratios estimated from galaxy col- 
ors but instead treat M*/L as an additional fitting parameter 
for each galaxy. 

There are two reasons why we do not adopt this "exact fit" 
method as our main approach. First, we choose the concen- 
tration parameter at random from the distribution P{c), and 
since the inferred disk mass fraction is sensitive to the value 
of c, the md values are not accurate on a galaxy-by-galaxy 
basis. Second, observational errors on M^,, Rd, and V2.2 in- 
troduce substantial uncertainty in the derived values of md, 
even if the assumed concentration is correct. (The same effect 
influences the A^ values in our standard approach, but much 
more weakly.) Despite these shortcomings, it is interesting to 
compare the P(md) estimated by this "exact fit" approa ch to 
the distribution assumed in the best-fit models of section l3^ 
Figure plots the inferred values of md against E*, as- 
suming AC and a Kroupa IMF (left panel) or a light IMF 
(right panel). Points show the md value determined for a halo 
concentration c = c(M). Error bars show the effect of raising 
or lowering the log of the concentration by iuc (0.14 dex). 
For higher c, the total halo mass must be lower to produce 
the same V2.2, and the value of md = M^,/M is consequently 
higher We have also tried varying the values of M*, Rd, and 
V2.2 by their 1 a observational errors, and these variations on 
their own produce an md spread similar to that shown by the 
error bars in Figure[T2l 

Lines in the left panel of Figure^jshow the mean trend and 
1 a scatter of our best-fit Kroupa-IMF model from ^3.21 The 
trend of the "exact fit" points is qualitatively similar to that of 
the model. However, there are a number of galaxies that can- 
not be fit with an md value below the universal baryon fraction 
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Fig. 12. — Disk mass fractions that reproduce the observed rotation velocities exactly for models with AC and two stellar IMFs: Kroupa IMF (left panel), light 
IMF (right panel). Vertical error bars illustrate the uncertainty associated with the scatter of halo concentrations. The end points correspond to logc ± av, with 
(7c = 0.14. The value of is not allowed to exceed 0.5. Lines show the mean trend (solid) and 1 a scatter (dashed) of the best-fit models with AC (eq. 1131 ^ for 
each IMF. 



of 0. 17, for any value of c in the expected range. (We have ar- 
bitrarily imposed an upper limit of 0.5 on md in Figure [T2l) 
These are the "maximal disk" galaxies from Figure^] which 
cannot be fit by an adiabatically contracted NFW halo with 
typical concentration. If we omit adiabatic contraction, some 
but not all of these points move to m,/ < 0.17. For the light 
IMF, only four galaxies require m,; > 0. 17 for a concentration 
in the expected range. The smaller number of physically pe- 
culiar results in this case is a circumstantial argument in favor 
of the light IMF 

It is tempting to go further and examine, for example, the 
correlation between 1114 and (see Jimenez et al. 2003). Un- 
fortunately, we consider correlations derived in this way to be 
unreliable. Adopting an incorrect value of c for an individual 
galaxy leads to an incorrect nid, and the inferred value of 
for that galaxy adjusts in a correlated way to reproduce the 
observed disk scale length R^. If we take an input model that 
has no correlation between m,/ and A^, we find that this ef- 
fect produces a strong, spurious correlation between the "ex- 
act fit" values of mj and A^. A similar spurious correlation 
arises from observational errors, as we illustrate in FigurefOl 
For this test, we create an artificial input sample by choos- 
ing the constant values of Aj = 0.05 and md = 0.08 for all 
observed galaxies and recalculating Rd and V2.2 accordingly, 
while keeping the observed M, and the mean halo concentra- 
tion c(M). Then we generate a realization of the new sample 
with added random "measurement" errors of M», Rd, V2.2 and 
perform the "exact fit", assuming a halo concentration drawn 
from P(c). The inferred values of nid and A^ scatter along 

1/2 2/3 
the lines of \d oc mj at high nid and Aj cx mj at low m^, 

displaying a spurious correlation. It might be possible to in- 
vestigate the correlation of ntd and \d by fitting resolved rota- 
tion curves and thus deriving the halo concentration for each 
galaxy, but even in this case one should be cautious of the 
impact of observational errors. 

We have also investigated the dependence of inferred val- 
ues of nid and \d on S,. The underlying question is: Do 
low surface brightness galaxies have low density because of 



low rrid (feedback or inefficient star formation) or because of 
high \d (high angular momentum)? We use the radius resid- 
ual, /S.Rd, relative to the mean Rd-M^, relation as our mea- 
sure of surface density variation at fixed M,. For our best-fit 
AC, Kjoupa-IMF model we find only a weak anti-correlation 
between ARd and nid, but we do find a significant positive 
correlation between ARd and A,/ (correlation coefficient is 
0.63). Thus in our best-fit model, the low-density galaxies 
tend to have both lower nid (a feature built into the model 
by the md-T,^, correlation) and higher Xd, but the Xd varia- 
tions dominate the spread of surface density at fixed M» . We 
are unable to confirm or reject this trend using our "exact fit" 
procedure, because any correlation is washed out by obser- 
vational uncertainties and by scatter of halo concentrations. 
Figure[0]shows that the inferred values of nid and A^ are scat- 
tered similarly for the high and low surface density galaxies. 
Even when we create another artificial input sample with an 
encoded anti-correlation of Xd and E*, the recovered values 
of Xd are completely uncorrected with S,. 

4. ADDITIONAL CONSTRAINTS 

So far, we have constrained the distribution of disk mass 
fractions using only the TF relation and its residuals, or the FP 
relation. In this section, we compare the best-fit models from 
^3.21 to additional observational constraints from the stellar 
mass function of galaxies and from the extended mass pro- 
files around galaxies probed by weak lensing and satellite dy- 
namics. For brevity, we will refer to the three models shown 
in Figures |8] and |9] all of which incorporate nid oc S'*", as the 
AC, no-AC, and light-IMF models. 

Given our best-fit P{nid\M^,,Rd), we can go from a cos- 
mological model prediction of the abundance of dark matter 
halos to a prediction of the number density of disk galax- 
ies as a function of stellar mass. To make such predictions 
we use the halo mass function calculated by the ellipsoidal 
excursion set method of Sheth & Tormen (2002,) . which ac- 
curately describes the results of cosmological A^-body sim- 
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Fig. 13. — Scatter in the infeiTed "exact-fit" values of mj and resulting 
from observational errors in M, , Rj, V2.2 and from scatter of halo concentra- 
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tions. The dotted lines show two power-law trends, A oc mj~ and A oc mj , 
respectively. True values were set to mj = 0.08 and A^ = 0.05 (at the in- 
tersection of the dotted lines). Solid squares and open circles correspond to 
galaxies with E, above and below lO'^MQkpc-^, respectively. 



ulations fe.g.. iSprineel & et alJ 120051) J We adopt cosmo- 
logical parameters of a ACDM model, il^ = 1 - = 0.24, 
= 0.041, 0-8 = 0.74, n = 1, //o = 70 km s"' Mpc-\ favored 
by the WMAP 3 - year data set and a varie ty of other cosmolog- 
ical observables dSoergel & et al.l l2006'). Note, however, that 
our adopted c(M) is calibrated on simulations with higher erg, 
a point we will return to shortly. We use the models defined by 
equation il3i to convert the halo masses M to stellar masses 
M,. First, we use the mean mass-radius relation (eq. 13) 
to express the average surface density as a function of stellar 
mass, cx M^^^, and convert the density-dependent model 
mrf(S]*) to the mass-dependent fhd(M^). This procedure yields 
similar but not identical functions to those given by equation 
J12t . Then we integrate equation to obtain the mean 
disk mass fraction for galaxies of mass M»: (mrf)(M») = 
/ mjP{mj\Mf:)d\ogmj, with an upper limit = 0.15. Fi- 
nally, we convert stellar mass to the mean halo mass, (M) = 
M,/ {nid) (M,), and shift the halo mass function to a galaxy 
mass function, dn„.J,M^,)/d\viMf, = dnii-iioi{M))/dlnM. 

Our calculation implicitly assumes that each dark matter 
halo hosts one and only one disk galaxy, so it does not in- 
clude the contribution of satellite galaxies to the galaxy mass 
function. However, theoretical models and observational es- 
timates from galaxy clustering and galaxy-galaxy lensing all 
indicate that the fraction of galaxies that are satellites in the 
mass range of our sample is low, between 10% and 20% 



' A study bv lJenldns et alJt2001D finds that the excursion set method gives 
the best fit to the W-body mass function when the halo mass is defined us- 
ing a constant overdensity 180 times the mean matter density. We have 
adopted a higher virial overdensity, 101 /f2m 421, which results in lower 
virial masses. In order to bring our halo masses to the same scale as the 
ISheth & Tormen >2002) mass function, we assume a halo outer density pro- 
file p(r) cx: r'^ and find the radius where the mean overdensity drops to 1 80. 
This radius is 33% larger than our adopted Rvir> and the enclosed mass is 15% 
higher than our Myi, , for a concentration of galaxy-sized halos c = 12. There- 
fore, we have shifted the .Sheth & Tormen 12002) mass function by a factor 
1.15 in halo mass. 



jMandelbaum et al.ll2005t IZehavi & et'al]l2005t laTeng et alJ 
l2005ir 

Figure [21 shows the theoretical halo mass function, the 
galaxy stellar mass functions predicted by the AC, no-AC, 
and light-IMF models , and the observational estimate of 
this mass function by I B ell et akl ([2003 ). derived from the 
SDSS and 2MASS surveys. We have converted this es- 
timate from Bell et al.'s "diet-Salpeter" IMF to a Kroupa 
IMF by multiplying their galaxy mass parameters by 10"** '^. 
The AC model substantially over-predicts the observed mass 
function at all masses, even without accounting for satel- 
lite galaxies. The high values of nid required to match 
the TF relation with a Kroupa IMF and AC puts galax- 
ies of a given M, into relatively low-mass, abundant ha- 
los, leading to a high number density. This tension be- 
tween the galaxy luminosity function and the Tully-Fisher 
relation is a long-standin g problem for CDM-bas ed models 
of galaxy formation (e.g., Kauffm ann et aLil993t ICole et alJ 
[1594; Somerville & Primack 1999). Here we are essentially 
confirming this problem, but removing some sources of un- 
certainty by using stellar masses in place of luminosities and 
by applying TF constraints from a well defined galaxy sam- 
ple. Of course, some of the halos in this range (Vyir = 80-300 
km s"') might host galaxies with significant bulges, which 
would not be in our TF sample - about 50% of the galaxies 
in Pizagno et al.'s (2006, in preparation) full sample satisfy 
the bulge-to-total < 0.1 criterion used here. However, if we 
lower the predicted mass function to account for the fraction 
of galaxies with substantial bulges, then we should also lower 
the observed mass function by the same factor, leaving the 
same g ap between them. Reconciling the AC model with the 
iBell et al. ( 2003 ) estimate for a Kroupa IMF requires, instead, 
that the halos that do not host disk-dominated galaxies have 
drastically lower values of m^, making them "dim" or "dark". 

The dotted curve in Figure [H] shows the prediction of the 
no-AC model, which lies below the Bell et al. (2003) estimate 
near the knee of the mass function but rises above it at low and 
high masses. We do not view the discrepancy above 10" M© 
as a serious problem, because a large fraction of galaxies in 
this mass range are early-type systems, which could plausi- 
bly have nid values lower by a factor ^ 2, and because mas- 
sive halos have a larger fraction of their baryon mass in satel- 
lite galaxies and intergalactic gas, leaving less for the central 
galaxy. The tension below ~ 2 x 10'" M© is more sig- 
nificant, since our sample extends to 2 x 10^ M0 and most 
galaxies in this mass range are late types, but it is arguably 
within the systematic uncertainties of the model and the data. 
We therefore consider the no-AC model to be an empirically 
acceptable way to reconcile the TF and mass function con- 
straints, albeit one that requires abandoning a seemingly well 
established element of galaxy formation physics. 

The dashed curve in Figure shows the model with adi- 
abatic contraction and a light stellar IMF. Rather than plot 
a second version of the "observed" mass function, we have 
added 0.15 dex back to the model values of M,, thus allow- 
ing a direct comparison to the original observational curve. 
If changing the IMF simply lowered nij values by 0.15 dex, 
then the predicted and observed mass functions would shift by 
the same factor, and the agreement between them would not 
improve. However, because the light IMF reduces the effect 
of disk gravity (both the direct contribution and the magni- 
tude of AC), it lowers the ratios V2.2/Kir and reduces average 
nid values by nearly a factor of two. This reduction improves 
the agreement with observations by putting galaxies in halos 
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Fig. 14. — Predicted mass functions of disk galaxies for the best-fit models 
with fhj oc EJ- with AC {dot-dashed line), without AC {dotted line), with AC 
and light stellar IMF {dashed line), and with low c{M) {long-dashed line). 
These galaxy number densities are calculated from the abundance of dark 
matter halos in the concordance ACDM cosmology as a function of galaxy 
stellar mass, M* = m^M, with the model values of mj{M) in each case. For 
comparison, top solid hne (labeled "halos") shows the case if the halos con- 
verted all of their baryons into stars. Bottom solid curve (labeled "galax- 
ies") shows the observed mass function of galaxies in the local universe from 
IBell etalJ 1200 3V The plotted range of masses con'esponds to the range of 
rotation velocities 80 — 300 km s"' , covered by our observed sample. 



of higher mass and lower abundance. The light-IMF model 
with AC lies midway between the Kroupa-IMF m odels with 
and without AC. It is marginally consistent with the lBell et alJ 
(|2003) mass function. 

A third way to reconcile the TF and mass function con- 
straints is to change cosmological parameters, and thus the 
halo mass function. We find that lowering the matter den- 
sity parameter to « 0. 15 brings the prediction of Kroupa- 
IMF, AC model low enough to match the knee of the observed 
mass function, and it also brings the light-IMF, AC model into 
rough alignment with the Kroupa-IMF, no AC model. Chang- 
ing the power spectrum normalization, o-g, on the other hand, 
has little impact on the predicted n(M,) because the halo mass 
function is only weakly sensitive to erg in this mass range. 

This discussion of the effects of changing cosmological 
parameters ignores the potentially crucial impact of ag on 
halo concentrations. Bullock etal. (2001b) present an an- 
alytic model for predicting halo concentrations (see also 
[Kuhlen et al. 2005), calibrated against their 57,^ = 0.3, erg = 1 .0 
ACDM simulation. If we apply this model to the WMAP 
3-year cosmological parameters, we obtain concentrations 
lower by about 25% at fixed mass, mainly because of the 
lower value of erg- In the range M = 10^' - lO'-'/i'^M©, the 
analytic model prediction is well fit by a power law, c{M) w 
9(M/1O'^/!"'M0)"° '. We have not adopted this as our stan- 
dard relation for the calculations above because the erg depen- 
dence of halo concentrations has not been thoroughly tested 
against A^-body simulations. However, the preliminary tests 
that have been carried out to date, and the tests of the pre- 
dicted redshift dependence of concentrations, suggest that the 
analytic model is accurate (J. Bullock, private communica- 
tion). 
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Fig. 15. — Average mass enclosed within 50, 100, 200, and 400 kpc (bot- 
tom to top) of galaxy centers as a function of stellar mass for the best-fit 
models with ifij oc Ej^-, with AC {solid lines), without AC {dotted lines), 
and with AC and fight stellar IMF {dashed lines). Pentagons show the aver- 
age halo virial mass for the standard IMF, with AC {filled) and without AC 
(open). 



The last lines of Tables ^ and |2] list TF and FP results for 
a model similar to our best-fit Kroupa-IMF, AC model with 
rhd cx S'J *'^, but using this lower c{M) relation. Lower con- 
centrations place a smaller fraction of the halo mass within 
2.2Rd, and therefore the total halo mass must be larger to 
reproduce the observed V2.2. The best- fit value of nid.Q (de- 
fined at M„Rf = lO^■^M0kpc"^ see eq. (l3l) drops from 
0.1 to 0.04. For a galaxy with the mean c(M) that lies on the 
mean Rj-M^, relation, this change of normalization keeps the 
mass of the adiabatically-contracted NFW halo within 2.2Rd 
approximately the same as before. Lowering c(M) and m^.o 
yields a slightly better fit to the TF relation, with weaker 
residual anti-correlation but slightly worse slope and inter- 
cept, and a noticeably better fit to the FP relation. However, 
the improvement in the predicted mass function is dramatic, 
as shown by the long-dashed curve in Figure^^ because the 
lower nid values put galaxies in more massive, less abundant 
halos. The low c(M), AC model fares just slightly worse than 
the high c(M), no-AC model in the mass function test. 

Using our best-fit models, we can also make predictions for 
extended mass distributions around galaxies like those in our 
sample. Figure^jshows the mean total mass enclosed within 
a fixed distance from the galaxy center, for distances rang- 
ing from d = 50 kpc to d = 400 kpc. For each galaxy in our 
sample, we use one realization of each model to calculate the 
interior mass Mint(< d) as a sum of the observed stellar mass 
and the predicted dark halo mass with the NFW profile. The 
distribution of gas (if present) is implicitly assumed to fol- 
low the distribution of dark matter At these large distances, 
we ignore the effect of AC on the shape of the NFW profile, 
although the model values are of course determined tak- 
ing AC into account at 2.2Rd. We average Mint(< d) in five 
bins of stellar mass containing approximately equal number 
of galaxies, and plot the mean Mint(< d) against the average 
stellar mass (M*) in the bin. For the AC and no-AC models, 
we also show the bin average of the virial mass (for critical 
overdensity A = 101), with the corresponding virial velocity 
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marked on the right axis. The interior masses for the no-AC 
model are higher than those of the AC model, typically by 
factors of 1.5-2, because of the lower values. The de- 
pendence of Mint(< d) on is shallower for the AC model 
because of the stronger dependence of on S,, which puts 
lower mass (and generally lower S,) galaxies in more mas- 
sive halos. The predicted interior mass in the models with 
fhd oc M*' is similar to that in the corresponding models with 
nid cx SC'. Predictions of the light-IMF model are close to 
those of the no-AC model, since they have similar values of 
niij. Predictions of the Kroupa-IMF model with AC and low 
c(M) are also similar to those of the plotted no-AC model at 
c/ < 100 kpc, but they lie systematically above it at larger d 
because the halos are more extended. 

These extended mass distributions can probed by galaxy- 
galaxy weak lensing observations, which yield the mean to- 
tal mass profiles after stacking a large number of galax- 
ies of similar luminosity or stellar mass. To compare with 
such measurements, we plot in Figure^] the virial mass-to- 
light ratios in r-band, Mvii /L,-, predicted by the AC models 
with the Kroupa IMF and the light IMF. For each galaxy in 
our sample, we draw values and concentrations from the 
P{md |M» , Rd) and P(c) distributions 1 times and compute the 
average Mvii, then divide by the galaxy's observed L,-. We 
follow the conventional practice of lensing studies in using 
Lr not corrected for internal extinction (which differs from 
standard practice in TF studies). Pentagonal points with er- 
ror bars show virial mass-to-light ratios of late-type galax- 
ies estimated by Mandelbaum et al. (2005), based on galaxy- 
galaxy lensing measurements in the SDSS. These ratios, rang- 
ing from ~ 40- 130Mq/Lq, are consistent with the light- 
IMF model but too high for the AC, Kroupa-IMF model. Star 
symbols with e rror bars sho w the galaxy-galaxy lensing re- 
sults of Hoekstr a et al J (120051) using the Red-Sequence Cluster 
Survey, which are somewhat lower and overlap with the pre- 
dictions of both models. Results for the Kroupa-IMF models 
with no AC and with AC but low c(M), not shown in the Fig- 
ure, are similar to those of the light-IMF model with AC. 

An independent measure of the dynamical mass at large dis- 
tances from the galaxy center is provided by the dynamics of 
satellite galaxies. Using a large sample of SDSS galaxies, 
Prada et al. (2006, in prep.) estimate Mvir = 2.4 x 10'- M© 
for the average luminosity = 1-8 x lO'^L© (open triangle 
in Figure [T6l. After converting luminosity to stellar mass, 
this corresponds to m,/ = 0.028, below our best-fit values at 
this M, {rhd = 0. 1 1 and 0.044, for the Kroupa-IMF and light- 
IMF models with AC, respectively, and = 0.042 for the 
no-AC model). Note that this estimated mass fraction ap- 
plies to an average of galaxies of all types, which could plau- 
sibly lie below the mass fraction of disk-dominated galax- 
ies in our sample. Stellar spheroids (bulges and elliptical 
galaxies) have higher stellar mass-to-light ratios than younger 
disk systems, so they have higher My^/L even for the same 
M,/Mvir. Mandelbaum et al. (2005j) find an average Myir/L,- 
of 140 for early-type galaxies vs. 53 for late types in the lumi- 
nosity range Lr = 10 - 10" L©. Furthermore, recent studies 
suggest that bulge-dominated galaxies have lower values of 
M^/Myir- For example. .Hump hrev et al. (2006) find a stellar 
fraction of 0.044 for an X-ray model of an elliptical galaxy 
NGC 720, while iHoekstra et al.. ilOOSH find on the average 
2.5 times lower M^/My„ for the red {B-V > 0.8) galaxies 
than the blue galaxies. Galaxies composed of a mixture of the 
disk and spheroidal components would have lower average 
M,/Mvir than pure disk-dominated galaxies, and this differ- 
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Fig. 16. — Dynamical mass-to-light ratio in the r-band, My\jLr, for the 
best-tit AC models with oc S'^- , for Kroupa IMF (filled squares and cir- 
cles) and light IMF {open squares and circles). In both cases, squares corre- 
spond to compact galaxies with i?j < RjiMt.) (eq. (3)), while circles con'e- 
spond to extended galaxies with Rj > Rj{Mt ). Plotted luminosities L,- are the 
observed values, not corrected for internal extinction. Star symbols with error 
bars show the galaxy-galaxy weak lensing measurements from Hoekstra et ^ 
1 2005), while open pentagons ai'e from Mandelba um et ali t200 5). Open tri- 
angle shows an average measurement using satellite dynamics from Prada et 
al. (2006, in prep.). 

ence may contribute to the gap between our AC models and 
the Prada et al. (2006) measurement. 

Available cosmological gasdynamics simulations of disk- 
dominated galaxies, which automatically include the halo 
response to baryon condensation, predict values in the 
same general range as our models. In a realiza tion of 19 
Milky Way-type galaxies, iSommer-LarsenI ( 12006 ) finds a 
mean baryon mass fraction = 0.08 at Vc = 200 km s"'. 
The TF relation for this simulated sample is roughly consis- 
tent with our observed relation ( Portinari & Sommer-Larse^ 
I2OO6I) . though the sample is too small to investigate residual 
correlations. For a single galaxy studied with ultra-high reso- 
lution, Robertson et al. (2004) a lso find nid » 0.08, whil e for 
a sample of three disk galaxies ^G overnato et al.l i2006h find 
values of mj between 0.04 and 0.08. 

All of our best-fit models introduce a correlation of mj with 
S*, which counteracts the tendency of disk gravity to pro- 
duce higher rotation speeds in more compact galaxies. For 
the model points in Figure ^] circles and squares respec- 
tively shows results forgalaxies above and below the mean 
Rd at a given (eq. |3l)- Larger R^ galaxies have system- 
atically higher M^t/Lr because they reside in more massive 
halos. If we divide the sample into three equal-size subsam- 
ples containing the smallest, intermediate, and largest galax- 
ies for their M*, the mean Kroupa-IMF My-„/Lr values are 
20, 23, and 28, respectively, and light-IMF Mvii /Lr values are 
42, 46, and 60, respectively. Galaxy-galaxy lensing analyses 
could test our proposed explanation of the weak correlation 
between TF residuals and radius, by detecting or ruling out 
this difference, roughly 40% between the smallest and largest 
1/3 of the population. 

Measurements of the rotation speed at distances larger 
than 2.2 Rd may also help to differentiate among our pro- 
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posed best-fit models. High density galaxies in our mod- 
els have relatively smaller dark matter halos and, therefore, 
more steeply falling rotation curves at large distances than 
low density galaxies. As an example, we consider hypo- 
thetical measurements, V3 2, of the rotation velocity at 3.2Rd. 
For the AC, Kroupa-IMF model, the average ratio V^^/Vi.i 
drops from 1 .00 for the lowest-density galaxies in our sample 
(M^Rf « 2 X lO^Mgkpc"^) to 0.97 for the highest-density 
galaxies (M ^ Rf Kil X lO^Mekpc"^). For the light-IMF and 
low c(M) models, V^^/yi.i varies from 1.02 to 0.99 in the 
same density range. In contrast, for the no-AC model the ra- 
tio V3.2/V2.2 is always in the range 1.00- 1.03, indicating that 
the rotation curve is still rising at 3.2Rd. Measurements at 
4Rd would show a larger difference. In the no-AC model, the 
ratio V4/V2.2 is still close to unity, while in the AC, Kroupa- 
IMF model it drops from 0.99 to 0.93, for the density range 
of our sample. This variation corresponds to up to 15 km s"' 
difference in the rotation speed, which could be detectable in 
a sample of extended, high-precision rotation curves. 

5. CONCLUSIONS 

We have constrained the structural properties of dark matter 
halos of disk-dominated galaxies using a well-defined sam- 
ple of 8 1 late-type galaxies with Ha rota tion curves, selected 
from the SDSS redshift survey dPiza gno et al. 2005). We 
model the Tully-Fisher and fundamental plane relations con- 
structed from the galaxy stellar mass, M», disk scale length, 
Rd, and optical Ha rotation velocity at 2.2Rd, ¥2.2- The stel- 
lar mass is determined from the /-band luminosity and a stel- 
lar mass-to-light ratio estimated from the galaxy's g-r color. 
The normalization of these stellar masses dep ends on the IMF , 
for which we generally adopt the form of 'Kroupal J2001I) . 
The maximum likelihood fit to the TF relation for our sam- 
ple is logy2 2 = 0.29(logM, - 10.5) + 2.23. The FP relation 
is log y2.2 = 0.27(logM,- 10.5) + 0.085(log/?rf- 0.65) + 2.23. 
Residuals from the TF relation are weakly, and positively, cor- 
related with residuals from the Rd-M^ relation. The weak- 
ness of this correlation, and the small R^ coefficient in the FP 
relation, imply that the TF relation is a nearly edge-on view 
of the fundamental plane of disk galaxies. 

We define models of the disk galaxy population by the prob- 
ability distribution function P(mij\M^,,Rd) of the disk stellar 
mass fraction nid = M,/Mvii-, where Myir is the virial mass of 
the host dark matter halo, assumed to have an NFW profile 
with the concentration predicted by A^-body studies. For each 
sample galaxy, we draw nid from the model distribution, then 
find the disk angular momentum parameter required to re- 
produce the galaxy's observed disk scale length. Our standard 
calculations include the impact of adiabatic contraction (AC) 
on the halo profile (e.g., Blumenthal et al. 1986; Gnedin et all 
^04), but we also consider models without AC. 

Models with a constant value of nid require nid ~ 0.08 with 
AC (rrid « 0.04 without AC) to reproduce the TF intercept 
and slope. These models under-predict the scatter in the TF 
relation, even when we include the full range of concentra- 
tion scatter found in A^-body simulations of dark matter halos. 
More importantly, they predict a strong anti-correlation of TF 
and Rd-M^, residuals, with more compact galaxies rotating 
faster at fixed . Equivalently, they predict an Rd coefficient 
in the FP relation that has the wrong sign and is incompati- 
ble with the observed value. Adding random scatter to V2,2 
values to match the observed intrinsic scatter does not resolve 
this conflict with the observed residual correlations. The con- 



flicts are more severe for the models with AC because of the 
stronger impact of disk gravity. The AC models also predict a 
TF slope that is slightly too steep. 

We conclude that a scatter of nid values, with dispersion 
in Innid of a,,, w 0.25, is required to reproduce the observed 
scatter of the TF relation, and that models with AC require 
higher disk mass fractions in higher mass galaxies to repro- 
duce the observed slope. These findings confirm the more 
qualitative arguments of Pizagno et al. (2005) based on the 
same data set. We can obtain consistency with the weak TF 
residual correlation and the small Rd FP coefficient by tying 
the mean disk mass fraction to the stellar surface density. The 
best-fit model with AC has iTid (x S"'^^ and a mean fraction 
iTid = 0. 1 at Rf = 10'^-^ Mq kpc"^ Without AC, the best-fit 
model has a shallower slope, 0.2 instead of 0.65, and a lower 
normalization, tfid = 0.04 at lO^'^Mo kpc~^. The m^-S, cor- 
relation counteracts the effect of stronger disk gravity in more 
compact galaxies. It produces a major improvement for the 
AC model, and a modest improvement for the no-AC model. 

The derived \d distributions for these two models are sig- 
nificantly narrower than the A distribution for dark matter ha- 
los in A^-body simulations: a\d ~ 0.25 (AC) or 0.39 (no AC) 
vs. a\ w 0.56. The mean Xd depends on the typical nid val- 
ues, which determine halo virial_ radii: we find Xd w 0.054 
(AC) or Xd w 0.036 (no AC) vs. A « 0.045 for A^-body halos. 
The systematic difference between the Xd and A distributions 
implies either that disk baryons have systematically different 
angular momentum from the dark matter in their host halo or 
that only a subset of halos form disk-dominated galaxies like 
those in our sample. 

We consider the m,/- correlation a physically plausible 
mechanism for bringing theoretically motivated disk forma- 
tion models, which include adiabatic contraction, into agree- 
ment with the slope and residual correlation of the observed 
TF relation. However, our best-fit AC model suffers the long- 
standing problem of over-producing the galaxy stellar mass 
function given the halo population of a ACDM cosmologi- 
cal model with observationally favored parameter values and 
the concentration-mass relation of Bullock et al. (2001b). It 
also predicts virial mass to r-band light ratios lower than some 
recent estimates from galaxy-galaxy lensing and satellite dy- 
namics. Lowering fim, and thus the halo mass function, can 
ameliorate the abundance problem, but the required value, 
r^m ~ 0.15, seems implausibly low and the mass-to-light ra- 
tio problem remains in any case. Eliminating AC can resolve 
both problems, at the expense of dropping a well-tested ingre- 
dient of galaxy formation theory. 

We suggest an alternative solution: adopting a "light" stel- 
lar IMF with M,/L lower than the Kroupa (2001) IMF by 
about 0.15 dex. This model yields acceptable residual corre- 
lations with a modest trend rhd oc E^J "^, predicts virial mass- 
to-light ratios in accord with recent measurements, and re- 
duces but does not entirely remove the tension with the galaxy 
baryon mass function. We do not know of other observa- 
tional evidence to support such a light IMF, but it is per- 
haps within the uncertainties of existi ng observ ational con- 
straints. Systematic uncertainties in the lBell et al.l (12003) pop- 
ulation synthesis models, which we use to compute M*/L 
from g-r color, could account for some of the required 
change even without altering the IMF itself. We also show 
that the lower c(M) relation predicted for = 0.74 by Bullock 
et al.'s (2001b) analytic model relieves most of the tension 
between the TF relation and the mass function and Myii/L,- 
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constraints, even in the case of a Kroupa IMF. A full assess- 
ment of this solution awaits better numerical confirmation of 
the analytic model and better observational constraints on 
Other changes to halo profiles that reduce the fraction of mass 
at small radii would have similar effects. 

Our observational findings are consistent with those of 
iDutton et al. (2006), which is reassuring given the different 
characteristics of the data sets. Our theoretical inferences 
are consistent but n ot identical. We are more reluctant than 
iDutton et a l.' ( 2006) to abandon adiabatic contraction, given 
the numerous simulations in which it occurs and the absence 
of any simulations in which it does not. The m^- corre- 
lation can remove the A7?^ - AV2.2 anti-correlation that other- 
wise plagues models with AC, though it leaves problems with 
the baryon mass function and weak lensing mass-to-light ra- 
tios. The addition of a light IMF or the adoption of lower 
halo concentrations largely resolve these difficulties as well, 
and we consider some combination of these effects to be a 
plausible resolution of the observational tensions. Modeling 
of rotation curve data at larger and smaller separations and 
galaxy-galaxy lensing measurements for more finely divided 
galaxy sub-classes may provide further insight that could help 



discriminate among these solutions. The fact that disk forma- 
tion models based on the most straightforward theoretical and 
observational inputs do not reproduce all aspects of the cur- 
rent data suggests that interesting surprises may still lie ahead. 
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APPENDIX 

MULTIDIMENSIONAL LINEAR FIT USING MAXIMUM LIKELIHOOD 

We have developed the following procedure for a multidimensional linear fit to data with observational errors in all variables, 
allowing for intrinsic Gaussian scatter. Each data point (galaxy) is represented by n + 1 parameters, «-vector x and scalar y, and 
their associated measurement errors, cti and ay. We label each galaxy by subscript ; and fit a linear relation of the form 



y=a{x-XQ) + b+N{0,a), 



(Al) 



where a is also an «-dimensional vector and xq is an arbitrary offset. The term A^(0,(t) accounts for normally-distributed scatter 
that is intrinsic to the studied relation and is not accounted for by measurement errors. In the case of the TF relation, such scatter 
could be due to the variations of the mass-to-light ratio among the galaxies, scatter in the properties of their dark matter halos, 
etc. We define the likelihood function 



2 In £ = - ^ In (cr^ + crj + a ,) - ^ 



[a{xi-XQ) + b-yi\^ 



+ const 



(A2) 



and find its maximum using Powell's direction set method (iPress et al.ll992h . The direction search is done as a function of « + 1 
variables (a and cr), since the value of the intercept b can then be calculated analytically: 

, Y.iWi[yi-a{xi-XQ)] 



Wi 



(A3) 



We estimate the errors of the best-fit parameters (Aa, Ab, Act) by bootstrapping the data. The error estimates are taken to be 
the standard deviation of the corresponding distributions of the best-fit parameters. We found empirically that of the order 1000 
bootstrap realizations are needed in order to obtain robust estimates of the errors, accurate to 3%. 

The offset xo can be chosen so as to minimize the uncertainty /S.b and eliminate correlation between Aa and A/?. The slope 
of the linear fit (a) and its error (Aa) do not depend on the choice of xq because its variation can be absorbed in the value of the 
intercept b. The intrinsic scatter and its error are also unaffected by the choice of xq. The error of the intercept b, however, does 
depend on xo and can be minimized by the appropriate choice of the offset. We use Brent's parabolic interpolation method (if 
« = 1 ) or downhill simplex method (if n > 1 ) to find the minimum of /S.b. 

Computer code that calculates the best-fit parameters and their error estimates is available from the authors upon r equest. An 
independent derivation of this maximum likelihood method and further discussion can be found in lD'Agostinilll2()()l . 
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TABLE 1 
TULLY-FlSHERFIT PARAMETERS 



Mi 


Pi 


Cm 




AC 




a b 




r 


X 


data 










0.291 ±0.013 2.231 ±0.006 


0.049 ±0.007 


0.23 ±0.14 




0.02 










N 


0.280 


2.265 


0.009 


-0.02 


67 


0.02 










Y 


0.292 


2.306 


0.014 


-0.20 


171 


0.04 










N 


0.295 


2.232 


0.016 


—0.26 


36 


0.04 










Y 


0.307 


2.275 


0.023 


—0.41 


84 


0.08 










N 


0.311 


2.196 


0.027 


-0.48 


68 


0.08 










Y 


0.322 


2.240 


0.034 


-0.57 


46 


0.04 







0.14 


N 


0.297 


2.232 


0.041 


-0.22 


13 


0.04 







0.14 


Y 


0.308 


2.276 


0.039 


-0.36 


70 


0.08 







0.14 


N 


0.313 


2.196 


0.040 


—0.43 


55 


o!o8 







0.14 


Y 


0.323 


2.240 


0'041 


-0.53 


40 


0.04 


+<T 





0.14 


N 


0.303 


2.233 


0.049 


-0.27 


14 


0.04 


+0- 





0.14 


Y 


0.314 


2.276 


0.049 


-0.42 


72 


0.08 
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0.14 
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0.14 
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-0.17 
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0.3 


0.21 


0.14 
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0.046 
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0.14 
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-0.07 
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0.65 
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0.041 


-0.18 


14 
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0.24 


0.14 


N 
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0.051 


+0.02 


6.9 


0.1 
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0.2 


0.14 


Y 


0.264 


2.247 


0.044 


-0.08 


16 














mrf oc E^^ Hght IMF 








0.04 


0.4 


0.25 


0.14 
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0.273 


2.275 


0.051 


+0.01 


4.6 














villi low c(M) 








0.04 


0.65 


0.2 


0.14 


Y 


0.279 


2.246 


0.041 


-0.06 


12 



Note - The intercept b is for the offset logM,o = 10.5. Lines with "+(t" correspond to models with added measurement errors to match the 
observed scatter In the ninth column, r is Pearson's linear correlation coefficient. We estimate its error by bootstrapping the observed sample. 

We find that Ar depends mainly on the number of objects in the sample, A^, and approximately equals W^^'^ . For the light IMF, the observed 
intercept h = 2.231+0.15 x 0.291 = 2.275. 
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TABLE 2 

Fundamental plane fit parameters 
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data 










0.271 ±0.017 0.085 ±0.054 
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0.04 










N 


0.312 


-0.075 


2.230 


0.008 


55 


0.04 










Y 


0.336 


-0.125 


2.272 


0.007 


111 


0.08 










N 


0.347 


-0.154 


2.192 


0.008 


125 


0.08 










Y 


0.369 


—0.200 


2.234 


0.007 


104 


0.04 







0.14 


N 


0.315 


-0.076 


2.230 


0.038 


18 


0.04 







0.14 


Y 


0.338 


-0.125 


2.272 


0.032 


77 


0.08 




Q 


0.14 


N 


0.349 


-0.155 


2.192 


0.030 


93 


0.08 







0.14 


Y 


0.370 


-0.201 


Z235 


0.024 


76 


0.04 


+(T 





0.14 


N 


0.325 


-0.093 


2.231 


0.048 


21 


0.04 


+(T 





0.14 


Y 


0.349 


-0.143 


2.273 


0.045 


80 


0.08 


+cr 





0.14 


N 




__n 17/1 


2.192 


0.043 


96 


0.08 


+(J 





0.14 


Y 


0.382 


-0.220 


2.235 


0.038 


76 


0.04 




0.22 


0.14 


N 


0.315 


-0.075 


2.230 


0.047 


16 


0.08 




0.27 


0.14 


Y 


0.361 


-0.172 


2.244 


0.038 


56 














fhd oc M£' 








0.04 


0.1 


0.25 


0.14 


N 


0.298 


-0.064 


2.234 


0.051 


10 


0.1 


0.3 


0.21 


0.14 


Y 


0.332 


-0.177 


2.235 


0.035 


41 














rhd oc 








0.04 


0.2 


0.25 


0.14 


N 


0.287 


-0.018 


2.234 


0.051 


4.6 


0.1 


0.65 


0.2 


0.14 


Y 


0.296 


-0.066 


2.241 


0.039 


14 


0.04 


0.4 


0.24 


0.14 


N 


0.261 


+0.022 


2.237 


0.051 


2.2 


0.1 


0.9 


0.2 


0.14 


Y 


0.270 


-0.022 


2.247 


0.043 


9.6 














mrf oc E^^ Hght IMF 








0.04 


0.4 


0.25 


0.14 


Y 


0.268 


+0.018 


2.275 


0.051 


1.8 














ihd oc S*^, low c(M) 








0.04 


0.65 


0.2 


0.14 


Y 


0.283 


-0.017 


2.246 


0.041 


9.8 



Note - The intercept b is for the offsets logM«o = 10.5 and logif^o = 0.65. For the Ught IMF, the observed intercept b = 2.275. 



Somerville, R. S. & Primack, J. R. 1999, MNRAS, 310, 1087 

Sommer-Larsen, J. 2006, ApJ, 644, LI 

Spergel, D. N. & et al. 2006, astro-ph/0603449 

Springel, V. & et al. 2005, Nature, 435, 629 

Steinmetz, M. & Navarro, J. F. 1999, ApJ, 513, 555 

Strauss, M. A. & et al. 2002, AJ, 124, 1810 

Strauss, M. A. & Willick, J. A. 1995, Phys. Rep., 261, 271 

Tully, R. B. & Fisher, J. R. 1977, A&A, 54, 661 

van den Bosch, F C. 2000, ApJ, 530, 177 

Vitvitska, M., Klypin, A. A., Kravtsov, A. V., Wechsler, R. H., Primack, J. R., 

& Bullock, J. S. 2002, ApJ, 581, 799 
Wechsler, R. H., Bullock, J. S., Primack, J. R., Kravtsov, A. V., & Dekel, A. 

2002, ApJ, 568, 52 



Weinberg, D. H., Colombi, S., Dave, R., & Katz, N. 2006, astro-ph/0604393 
White, S. D. M. & Rees, M. J. 1978, MNRAS, 183, 341 
York, D. G. & et al. 2000, AJ, 120, 1579 

Zavala, J., Avila-Reese, V., HemSndez-Toledo, H., & Firmani, C. 2003, A&A, 
412, 633 

Zehavi, I. & et al. 2005, ApJ, 630, 1 

Zheng, Z., BerUnd, A. A., Weinberg, D. H., Benson, A. J., Baugh, C. M., 
Cole, S., Dav^, R., Frenk, C. S., Katz, N., & Lacey, C. G. 2005, ApJ, 633, 
791 



